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ABSTRACT 

Stars with helium cores between ~64 and 133 Mq are theoretically predicted to die 
as pair-instability supernovae. This requires very massive progenitors, which are the- 
oretically prohibited for Pop II/I stars within the Galactic stellar mass limit due to 
mass loss via line-driven winds. However, the runaway collision of stars in a dense, 
young star cluster could create a merged star with sufficient mass to end its life as a 
pair-instability supernova, even with enhanced mass loss at non-zero metallicity. We 
show that the predicted rate from this mechanism is consistent with the inferred vol- 
umetric rate of roughly 2 x 10~^ Mpc^"^ yr^^ of the two observed pair-instability 
supernovae, SN 2007bi and PTF lOnmn, neither of which have metal-free host galax- 
ies. Contrary to prior literature, only pair-instability supernovae at low redshifts z < 2 
will be observable with the Large Synoptic Survey Telescope (LSST). We estimate the 
telescope will observe 10^ such events per year that originate from the collisional 
runaway mergers in clusters. 
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1 INTRODUCTION 

Pair-instability supernovae (PISNe) are thought to occur 
for stars wit h helium cores between ~64 and 133 Mq 
l|Heger fc Woo slcv 200^). At zero metallicity, this corre- 
sponds to initial stellar masses between ~ 140 and 260 
Mq . These enormous stellar masses may have been reached 
by Pop II I stars, predicted to ha ve a top-heavy mass dis- 
tribution l|Bromm fc Larson|[2003 '). However, at lower red- 
shifts, as the universe was enriched. Pop III stars ceased to 
form once the local metallicity exceeded a critical thresh- 
old Zcrit ~ W'^Zq (|Bromm fc Loebl 120031 '). Since it is al- 
most impossible to raise the int ergalactic medium metal- 
licity in a homo geneous way (|Furlanetto fc Loebl l2003l : 
IScannapieco et al]i2003 ). pristine metal- free stars will still 
be formed past th e end of the reionization epoch z < 6 
Trenti et al.ll2009l). co nceivably all the way down to z = 2.5 



Tornatore et al.lboOTi ). Observations have confirmed the ex- 



istence of extremely metal-poor star fo rmation at moderate 
redshifts of z = 3.357 and z = 5.563 (|Fosburv et al.ll2003l : 
iRaiter et ailbOld '). The detectability of PISNe from Pop 
III stars at these mode rate redshifts was investigated by 
IScannapieco et al] (|2005l '). 

Outside these surviving pristine regions, there is a wide 
range of observations that support an upper limit to stellar 
mass at ~ I5OM0 in our Galactic neighborhood (iFigeijii2005|; 



IWeidner et al.l [20101 ) . prev enting the formation of PISNe 
from Pop II/I stars (but see lCrowther et al] (|2010l ) for stars 
determined to be above I5OM0 in the R136 cluster.) Nev- 
ertheless, even if very massive stars can form in metal rich 
regions, these radiatively supported stars are loosely bound 
and have strong winds driven mainly by radiati on pressure 
through spectral lines , scaling as M oc ^''.5~o.7 (jvi^i; gt al.l 
l200ll : lKudritzkill2002h . So even Pop II/I stars with initial 
masses between ~ 140 and 260 Mq will suffer copious mass 
loss during both the hydrogen and helium burning stages, 
and may not end their lives with enough mass remaining to 
die as PISNe; this prediction could be contested, as there 
are still large uncertainties in mass loss models from hot 
massive stars (|Puls et al.ll200^ ). Nevertheless, due to mass 
loss the possibility of PISNe is usually not considered for so- 
lar composition stars, even though the pair instability arises 
irrespective of the progenitor's metallicity. 

Regardless, PISNe have very likely already been ob- 
served at low redshifts, most convincingly in the case of 
the very luminous and long duration event SN 2007bi 
()Gal-Yam et al.ll2009l ) . More recently, the Palomar Transient 
Factory observed a new presumed PISN, PTF lOnmn (Gal- 
Yam 2011, submitted; Yaron et al., in preparation), and an- 
other PISN candidate was reported by the Pan-STARRSl 
survey, PSl-llap (Kotak et al., in preparation.) Although it 
may be possible to explain bright events like SN 2007bi with 
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altern a tive models (e.g. Woosley et al.l (|2007l ) , iMoriva et al.l 
l|2010l '). lKasen fc BildsterJ (|201(]| ')') on the whole the observa- 
tions seem to favor a scenario in which a large total mass 
and radioactive mass were ejected, as in a PISN explosion. 
The observations therefore suggest that very massive stars 
above the Galactic limit are formed in the local universe. 
The metallicities of the host galaxies of both supernovae 
are low but well abo ve the maximum m etallicity required 
to form Pop III stars (| Young et al.ll2O10l ). Either pockets of 
pristine gas survived in the dwarf host galax;y of SN 2007bi 
and FTP lOnmn, or the initial mass function (IMF) of Fop 
II/I stars merely steepens at the very high end (instead of 
a hard upper limit), or there are other exotic ways to form 
a very massive star. 

In theory, mergers of stars can form massive SN progen- 
itors at any metallicity and circumvent the upper mass limit 
for Fop II/I stars at ~ ISOM©. The most likely environment 
for such mergers is a dense, young star cluster undergoing 
core-collapse, in which a runaway collision product can be- 
come massive eno ugh to die as an ultra-luminous supernova. 
iForte gics Zwart fc van den H cuvcl (2007) first inve s tigate d 
this scenario fo r a co U apsar, and lYungelson et al. ] (|2008l ): 
iGlebbeek et all (|2009h : IVanbeveren et al.l (|2009l ') discussed 
the conditions under which the runaway collision product 
will end its life as a FISN. 

In this paper, we calculate the number of the collision 
runaway merger products within dense young star clusters 
that lie in the FISN progenitor mass range, and show that 
the predicted event rate is roughly equal to the inferred rate 
of FISNe from the detection of SN 2007bi and FTP lOnmn in 
existing surveys, without requiring the supernova progenitor 
to be metal free. We further investigate the observability 
and rate of these events in the low redshift universe with 
the Large Synoptic Survey Telescope (LSST). 



2 RATES FROM RUNAWAY COLLISIONS 



An ap preciable fraction of stars are born in clusters; iBastianI 
l|2008l ) found the fraction of mass that forms in clusters > 
IOOMq out of the total star formation rate to be P 8 ±3%. 
As soon as the cluster forms, the massive stars start to sink 
to the cluster center due to dynamic friction, driving the 
cluster into a state of core collapse on a timescale of tec ~ 
0.2trh, where t,.h is the relaxation time: 



trh 



»--(t^)'(t«% 
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Here m is the cluster mass, r is its half mass radius, TV is the 
number of stars, (m) = N/m ~ 0.5M p, is the average stellar 
mass , and log A « log(0.1A'') ~ 10 (|Fortegies Zwart et al.l 
I2OIOI ). In sufficiently compact clusters, the formation of a 
dense central subsystem of massive stars may lead to a 
collision runaway, where multiple stellar me rgers result in 
the formation of an unus ually massive object (|Giirkan et al.l 
l2004l : lFreitag et ahllioO^ ). This prescription is often invoked 
to form intermediate-mass black holes via the photodisin- 
tegration instability that collapses a super-massive star di- 
rectly into a black hole. 
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Figure 1. Differential number distribution of tiie final runaway 
mass formed, per IMq of stellar mass formed in all clusters. The 
calculated distribution is not perfectly smooth owing to the finite 
number of samples in the observed radius distribution. 



Por a successful collision runaway to occur, the star 
cluster must experience core collapse before the most mas- 
sive stars explode as a SN (~3 Myr). Por compact clus- 
ters {trh ^ 100 Myr), basically all massive stars sink to 
the cluster core during the runaway, and the final merged 
object's mas s scales with the cluster mass, ~ 8 x 
10"^mlogA jPortegies Zwart fc McMillan|[2003 ). Por clus- 
ters with longer relaxation times, only a portion of mas- 
sive stars sink to the core in time and the merged object's 
mass scales as mt~^^^ ( ,McM illan & Fortegics Zwart 200^1). 
A fitting formula for combi ning these scalings is given by 
iFortegies Zwart et al.l l|2006l ), calibrated by N-body simula- 
tions for Salpeter-like mass functions: 



rrir ~ 0.01m(l -|- 



-)" 



(2) 



lOOMyr ' 

To get statistics on the final runaway mass rrir from equa- 
tions ([T]) & ((21), we need to specify the number distribution 
of clusters as a function of their mass m and radius r. 

The functional form of the clu ster in itial mass function 
is well represented by a lSchechtei] (| 19761 ) distribution. 
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~ 2 (|Zhang fc Palll Il999l : 
Por Milky Way-type spira l 
galaxies the break mass m« « 2 x 10^ l|Gieles et al II2OO6I : 
lLarsen|[2009l ) , whereas for i nteracting gala, xies and luminous 
IR galaxies m* > WMq (iBastianI (200I) . Our resuhs are 
not sensitive to the choice of m*. 

Several studies have discussed the lack of any clear cor- 
relation between the s i ze of a cluster and its m ass or lu- 
minosity l|Larsenl 12004 IScheepmaker et al.l |2007| ). Lacking 
a functional distribution of cluster radii, we use the em- 
pirical distribution of radii for each cluster mass bin, for 
obser ved clusters younger than 5 Myr, compiled in Tables 
2-4 of lPortegies Zwart et all (|20ld '). The restriction on clus- 
ter age is important, as clusters expand considerably during 
the first 10 Myr of their evolution. Note that this empirical 
construction underestimates the number of super-massive 
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Figure 2. Predicted rate of PISN events per comoving Mpc'^ 
per year. The pair-instability SNe progenitor mass range is a 
major uncertainty for non-pristine stars. Here we use 140-260 
Mp, (blue line) for meta l-free Pop III stars from the models of 
iHeger fc WooslevI whe reas the mass range of ~250-800 

Mq (red line) is taken from lYuneelson et al.l [200Sl ). who ac- 
count for increased mass loss at solar metallicity. The environ- 
ments of low redshift PISNe will likely lie in betw een these two 
cases. A stronger mass loss scenario is presented bv lBelkus et alj 
1I2OO7I ). who found that PISNe progenitors can only be created at 
metallicities below 0.02 Zq, with a mass range ~300-1000 Mq 
(black line); as the fraction of matter in Z < 0.02 Zq galaxies 
is not well constrained past low rcdshifts, wc do not plot this 
rat e past z = 0.5. The strongest mass loss scen a rios p resented 
by iGlebbeek et al.l (2009) and Vanbeveren et al.l l l2009b predict 
a PISNe rate of practically zero, so we do not plot i t here . The 
dashed lines are the PISNe rates of lYungelson et al.l l|2008h and 
iBelkus et al. I ||2007|) adjusted for mass loss from stellar collisions. 
The black star shows the inferred PISN rate from current surveys. 



collision runaway objects (> 10^ Mq), as there happens to 
be no observed > W^Mq clusters younger than 5 Myr in 
the current sample, but this does not drastically affect our 
PISN rate estimates. With the joint number distribution of 
clusters as a function of their mass and radii, we can find 
the number distribution of the final mass of the runaway 
collision merged object, see Figure [1] 

However, as we have not taken mass loss into account 
in our estimate of the final runaway mass in equation 
we artificially inflate the mass range for PISN progenitors 
required at the end of the last merger event, to compen- 
sate for the mass lost during the collision runaway merger 
sequence. For zero-metallicity Pop III st ars, mass loss via 
line- d riven winds should be negligible, and lHeger fc WooslevI 
l|2002ri found the progenitor mass range to be 140-260 Mq; 
this should set the upper limit of the PISNe rate from run- 
away collision products. 

As for Pop II/I PISNe progenitors, we caution that mass 
loss via stellar winds for massive stars M > IOOM0 is not 
well underst ood. In fact, the o bservations of PISNe at low 
redshifts (G al- Yam et al.ll2009l ). and of Type Iln SNe whose 
progenitors are found to sometimes retain th eir hydrogen 
envel opes until shortly before their explosion (|Smith et al.l 
I2OIII ). suggest that most commonly-used stellar mass loss 
models are inaccurate for very massive stars, and likely over- 
estimate the total mass loss, as the models do not allow such 



Figure 3. Number of new PISNe per deg^ per unit redshift per 
year, for the same models as Figure[2] Note that LSST is expected 
to cover over 20,000 deg^ of sky. 



SNe to exist at the measured metallicities. Therefore, to ac- 
count for this uncertainty, we present here various PISN 
progenitor mass range scenarios described in literature, de- 
pendent on the assumed metallicity and mass loss prescrip- 

tioU; 

lYungelson et al. 1 (|2008l ) studied the evolution and 
fate of super-massive stars with solar metallicity from the 
zero-age main sequence using detailed stellar structure 
models. However, instea d of extrapolat i ng co m monly used 
mass lo ss models, e.g. de Jager et al.l (Il988l ): Vink et all 



l|200ll) : iKudritzkil (|2002l ). lYungelson et al.l (|2008l ) used 



an ad-hoc mass-loss prescription consistent with existing 
models in their relevant regimes and more consistent 
with the observed Hertzsprung-Russell diagram location 
and mass loss ranges found for young massive stars in 
clusters in the Milky Way and the Magellanic Clouds. 
Notably, their time-averaged Wolf-Rayet (WR) mass loss 
rate Mwr hardly exceeds 10"* Mq yr~^, which better fits 
observations of hydrogen-rich WR stars that account for 
i ron-line blanketing and clumping in determining Mwr 
( Ha mann et al.l l2006 h. and also agre es well with Mw r 
estimates based on radio observations (|Cappa et al.|[20o3 ). 
On the contrary, the extrapolation of Wolf- Rayet mass los s 



formulas to high st ellar mas ses given by Lange i^ (Il989|), 
iNug is fc Lamerd (120001). [Nelemans &: van d en Heuvel 



VanbeverenI ||2003l ) overestimate 



(2001'), and iDe Ponder 
the mass loss rates compared with these observations . 
Therefore, we use the results of lYungelson et all (I2OO8I ) 
as our fiducial model. They allow the creation of PISNe 
progenitors at Z ~ Zq in the initial mass range of ~250-800 
Mq; however, they do not account for the mass loss from 
stellar collisions. 

Alternatively, by extrapol ating theoreti c al m ass-loss 
rates for radiation-driven wind, iBelkus et al.l (|2007l ) found 
that when the metallicity Z is between 0.001 and 0.02 Zq, 
one may expect PISN candidates for stars with masses from 
~300-1000 Mq; however, at Z > 0.02 Zq no PISNe are 
expected. Usi ng the observed ga l axy luminosity-meta llicity 
relationship (|Kirbv et al.l l2008l : iGuseva et all |2009| ) and 
the galaxy luminosity function at low redshifts ^ ~ 0.1 
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l|Blanton et al.ll2003h . we find that ~ 0.3% of stellar mass 
is formed in Z < 0.02 Zq galaxies, and fold this factor into 
the predicted PI SNe rate for this scenario. 

In addition, iGlebbeek et"al] l|2009l ) follows the evolu- 
tion of the collision product for a few merger sequences 
for a m ~ 5 X 10^ A/© cluster, including mass loss along 
the co urse of the c ollisio n sequence by using the prescrip- 
tion of IVink et all l|200ll 'l. and found that above Z = 0.001 
Zq, the collision runaway product cannot die with suffi- 
cient mass to undergo a PISN. The main sequenc e stellar 
wind mass loss rate betwee n this work and Yungcls on et al.l 
l|2008l ) are similar, however, iGlebbeek eral.l ((20091) also cal- 
culates the mass loss from stellar collisions to be roughly 
~ 20% of the total merger product mass before mass loss. 
Nevertheless, the main source of discrepancy between their 
conclusio ns is due to their very different Wolf- Ray et mass 
loss rates. IGlebbeek et al.l (|2009|) implements a stron g Wolf- 
Rayet mass loss rate from iNugis fc Earners! (|2000l ) (up to 
3.6 X 10"^ Mq yr"^ at Z = 0.02), bringing the colli- 
sion product down to only rrir ~ IOMq by the end of 
core helium-burn i ng. U sing a comparable mass loss rate, 
IVanbeveren et all (|2009l ) reaches the same conclusion that 
PISNe cannot occur above Z = 0.001 Zq. Note that with 
these mass loss rates, essentially no star in the low redshift 
universe below M ~1000 A^q will end their lives as a PISN, 
irrespective of the collision runaway mechanism. 

To account for the ~ 20% mass loss due to unbound 
ejecta from the stellar collision, we can further increase the 
required PISN progenitor mass range. The new adjusted 
mass range for PISN p rogenitors would be ~313-1000 Mq 



Yungelson et al. (|2008.) scenario, and ~375-1250 Mq 



Belkus et al.l (|2007l ) scenario 



in the 
in the 

Combining the above, we can estimate the number of 
collision runaway products that have a final mass rrir in the 
various PISN progenitor mass range sce narios. Using the 
lobal comoving star formation rate from iReddv fc Steidell 
20091 ). we estimate the PISN rate as a function of redshift 
in Figures [2] and [3] If the collision runaway mechanism is 
indeed responsible for creating PISNe progenitors at Pop 
II /I environments in the local univer se, we find that only the 
mass loss prescription described bv I Yungelson et al.1 (|2008l ') 
fits the current rate of PISNe inferred from observation. 



3 OBSERVABILITY WITH LSST 

The Large Synoptic Survey Telescope is a planned wide-field 
survey telescope that should begin operations at the end of 
this decade. It has a very wide field of view of 9.6 deg^, 
and 6 bands: u,g,r,i,z, and y, covering 320-1080 nm. For 
the most sensitive bands g, r, and i, a single visit will reach 
Mab = 25.0, 24.7, and 24.0 (5(t) sensitivity, respectively. 
These bands will be visited 10, 23, and 23 times every year 
during the 10 years of operation, reaching a coadded depth 
of Mab =26.3, 26.4, and 25.7 per year by stacking multi- 
ple images. Note that for objects much dimmer than ~ 22 
mag/arcsec^, or ~ 25.5 mag/pixel for LSST, the signal will 
be dominated by the sky background (e.g. airglow and zodi- 
acal light), so in this regime the limiting signal fiux needed 
to reach a fixed signal-to-noise ratio is inversely proportional 
to the square root of the integration time. 

We use simulated PISN light curves and spectra from 
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Figure 4. Maximum redshift observable by LSST, as a function 
of the number of stacked images, for various PISN progenitor 
models. Here we use the co-added r band 5cr sensitivities, for 
which LSST will visit the same location 23 times every year, or 
once every ~ 16 days on average. We consider a PISN at a certain 
redshift as observable if it stays brighter than the limiting co- 
added depth for a duration longer than the time it takes to observe 
that number of images. Zmax eventually drops with increasing 
image count, as the PISN flux falls off but the sky background 
remains, reducing the integrated signal-to-noise. The brightest 
PISN will be observable with LSST out to a redshift of ~ 1.8. 



iKasen et al.l l|201ll ). who improved radiative transfer calcu- 
lations by using a multi-wavelength Monte Carlo code which 
includes detailed line opacities. In particular, we use models 
Hel30, HelOO, and He80, which represent pair-instability ex- 
plosions of non-rotating bare helium cores with masses 130, 
100, and 80 Mq, respectively, as non-pristine massive stars 
formed via runaway collision will likely lose most of their 
hydrogen envelope by the end of their life. The brightest 
helium core model Hel30 peaks at around Mab ~ —22 in 
the rest frame r band, and stays above Mab ~ —21 for half 
a year, and above Mab — —20 for almost one year. Such 
an event in the local universe will be easily detectable; how- 
ever, the rates for PISNe from both Pop III and Pop II/I 
progenitors is predicted to be very low at 2; ~ 0. These rates 
increase at higher redshifts, but since the higher wavelength 
z, y LSST bands are much less sensitive, the best strategy 
to find PISNe is to continue using the g, r, and i bands and 
observe at the rest frame UV and optical luminosity of the 
supernovae. 

Using the coadded depth sensitivities, we find that us- 
ing the r band is optimal for the helium core PISN models, 
and that we can observe the brightest He 130 model out to 
a redshift of 2 ~ 1.8 by stacking ~ 10 images (see Figure 
|4|. Below z < 2, the PISN is visible in the r band for over 
1 year in the observer frame; however, at 2 ^ 2, the su- 
pernova will be too dim in the rest frame UV wavelengths 
being effectively probed, even though the (1 -I- 2) time di- 
lation allows more stacked images. Even if one combines 
data from the g,r, and i bands over one year, and reaches 
a coadded depth of Mab ~ 27, the supernova will still be 
too dim to be observable beyond 2 = 3. Alternative PISN 
models where the progenitors are red supergiants which re- 
tain their hydrogen envelopes have a longer plateau in their 
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light curves, and thus stay visible slightly longer than the 
helium core models. However, the conclusions are similar - 
in terms of instrument capability, redshifts z < 2 axe most 
suitable for detecting PISNe in the normal operation mode 
of LSST. The smaller HelOO model is only visible out to 
z ~ 1.2, while even smaller progenitors are too dim to be 
seen beyond z < 0.4. Combined with Figure [S] we estimate 
that LSST will see on the order of ~ 10^ new PISNe per 
year that originated from the final collision runaway object 

in young, dense clusters. 

These conclusions differ from those of iTrenti et al.l 

20091) as well as the L SST Science Book 

LSST Science Collaborations et all |2009| ). which con- 
cluded that PISN at jz ~ 4 will be w ithin the capability 
of LSST. The difference arises because iTrenti et all \20m ) 
approximated the PISN with a blackbody spectrum with 
Tsff — 1.5 X lO^K, which overestimates the rest frame UV 
flux compared to the spectru m obtained by the radiation 
hydrodynamics simulations of iKasen et al. 1 (|201lD . Also, in 
the LSST Science Book, when calculating that hundreds 
of z = 2 - 4 PISNe will be detected by LSST (Chapter 
11.14), the authors used z and y band sensitivies of ~ 26.2. 
This is unrealistic as Mab ~ 26.2 can only be reached in 
the z band by stacking all images over the entire 10 year 
lifetime of the survey, but no PISN will stay bright enough 
that long even with time dilation; the y band is even less 
sensitive. Our findings suggest that, to find PISNe at z > 2, 
an instrument with better infrared capabilities such as the 
James Webb Space Telescop^ is required. 

Although stacking multiple images averages out the 
time variation in the supernova light curve, LSST also al- 
lows a secondary survey over a smaller area of sky, going 
substantially deeper in a single epoch. However, due to the 
steep luminosity function of PISNe, we will preferentially 
see only the massive PISN events beyond the local universe, 
so narrow, deep exposures by LSST are more useful for im- 
proving light curve coverage, instead of supernova discovery. 



4 DISCUSSION 

Runaway collisions were explored most seriously in mas- 
sive, dense clusters, so equation ((2]) may not be accurate for 
m < 10'* Mq. However, only more massive clusters can make 
runaway masses rrir in the PISN progenitor mass range, so 
this does not affect the predicted PISN rate. In addition, 
initial m ass segregation of sta rs wi thin young c l uster s ob- 
served bv lde Griis et al.l l|2002h and lStolte erall (|2006l ) will 
shorten the time to runaway collisions and increase rrir, but 
we do not take this into account. 

For z < 6, the rate and detectability of PISN from 
Pop HI progenitors bo rn in surviving pockets of metal-free 
gas was investigated bv lScannapieco et al.l l|2005l '). To model 
the PISN light curves, they used an implicit hydrodynamics 
code which only implements gray diffusive radiation trans- 
port; for spectra and colors they assumed a blackbody dis- 
tribution. Depending on the intergalactic medium metal en- 
richment history, their predicted rates span two orders of 



http://www.jwst.nasa.gov/ 
© 2011 RAS, MNRAS 000, [iHT] 



magnitude, with their lower end roughly equal to our colli- 
sion runaway rates at z = 1 — 2. However, a PISN with a 
pristine host galaxy has yet to be observed. 

A pilot search done using the Spitzer /IRAC dark 
field found no candidates above the sensitivity limit of 
AfAB(3.6^m) ~ 24, placing an 95% confidence upper limit 
of 23 per deg^ per year for > l ^tjy sources with plateau 
timescales less than 400/(1 -I- z) (|Frost et al.l 120091 '). which 
does not contradict the predicted rate of < 0.1 PISN per 
deg^ per year for our collision runaway model. 

More recently, observers have discovered a class of 
ultra-luminous supernova, with luminosities exceeding those 
of the brightest pair-instability events, and rates of order 
~ 10~* Mpc~^ yr~'^at z ~ 0.3. These events do not ap- 
pear be standard radioactively powered PISNe, as their lu- 
minosities are too high and their light curve durati ons too 
short (e.g.. lQuimbv et al.|[2oTll : IChomiuk et aLlbOllI ). Com- 
paring the rate of those events to that of the two putative 
observed PISNe, G al- Yam found that PI SNe are roughly ~5 
times rarer than the lOuimbv et al] (|201lh ultra-luminous su- 
pernovae (Gal- Yam 2011, Science, submitted). This gives a 
PISN rate of ~ 2 x 10"'' Mpc"^ yr"* in the local universe, 
roughly consistent with the collision runaway rates found in 
Figured 

If the collision runaway of massive stars in young, mas- 
sive stellar clusters do give rise to PISNe at Pop II/I metal- 
licities, we expect to see such a young, massive cluster at 
the same location, after the light of the supernova has faded 
away. However, even a 10^ M© cluster only has an absolute 
magnitude of about -8.2 mag, so the PISN will have to oc- 
cur close by {z < 0.05) for its host cluster to be observed 
with current telescopes. Also, these PISNe should follow the 
distribution of clusters, and appear in the luminous parts of 
their host galaxies, analogous to the position of long dura- 
tion gamma-ray bursts. Note that due to the steep distri- 
bution of collision runaway masses (see Figure [TJ, the rates 
of PISNe from collision runaways will still be higher in en- 
vironments with low metallicities, as long as mass loss for 
massive stars is proportional to metallicity. 

Alternatively, if mass loss models are wrong and the 
Galactic stellar mass limit is violated, we need not invoke 
stellar mergers to create the massive progenitors r e quired 
for the observed non-pristine PISNe. iLanger et al. ] (|2007h 
found that hydrogen-rich PISNe could occur at metallicities 
as high as Zq/3, resulting in a rate of about 1 PISN per 
10'^ SNe in the z « universe. For a more conservative 
metallicity threshold of Zq/10, the rate would be about 1 
PISN per 10* SNe. However, even the latter is a few times 
higher than the current inferred rate of PISNe. 



5 CONCLUSION 

We have shown that the runaway collision and merger of 
stars in a young, dense star cluster may form the massive 
progenitor of a pair-instability supernova at non-zero metal- 
licity. The volumetric rate of such events is a few times 10~® 
Mpc~^ yr~^ in the local universe, roughly matching the in- 
ferred rate of pair-instability supernova events SN 2007bi 
and PTF lOnmn in ongoing surveys, both of which have a 
metal-poor but not metal-free host galaxy. We expect that 
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the primary survey of the Large Synoptic Survey Telescope 
would see ~ 10^ such events per year. 
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